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We apply the target rest frame formalism to photon + heavy quark production cross section in 
hadronic collisions at high energies. We investigate the dependence of the production cross section 
on the photon and quark rapidities and transverse momenta. It is shown that the photon transverse 
momentum spectrum is a sensitive probe of color dipole scattering amplitude. The theoretical results 
are compared to Tevatron measurements of the differential 7 + c+ X and 7 + 6+ X production cross 
sections at ^/s = 1.96 TeV. An analysis for proton-proton and proton-lead collisions at the LHC 
regime is also performed. 



PACS numbers: 12.38.Aw, 13.60.Hb, 13.85.Qk 

I. INTRODUCTION 

The inclusive direct (prompt) photon production in 
hadronic collisions has been quite useful for accessing 
information on parton distributions in hadrons. More- 
over, direct photons appears to be an useful probe of 
the initial state of matter created in heavy ion collisions. 
Their interaction with the medium is purely electromag- 
netic, which gives a baseline for the interpretation of jet- 
quenching models. From the theoretical point of view, 
the prompt photon production is reasonably described 
in the next-to-leading order (NLO) perturbative QCD 
(pQCD) approach fl| and on the color dipole formalism 
I^-lZi] as well. Recently, the first measurements of inclu- 
sive photon production in association with heavy flavour 
jets have been reported by DO Collaboration (8|. Such 
a process provide important information on the parton 
content of the initial state hadrons as it is sensitive to 
charm (bottom) and gluon densities within the colliding 
hadrons. 

In pQCD the inclusive production of prompt photons 
plus a heavy flavour quark 9] is driven by the QCD 
Compton scattering, g + Q — > 7 + Q, and also receives 
contribution from quark- antiquark annihilation process, 
7-l-(7 — 7-f QQ. This is the reason of such a pro- 
cess to be sensitive to the heavy quark and gluon content 
of hadron, which have considerable uncertainties. The 
data description was found to be fairly good [1] for large 
transverse momenta photons, where a perturbative ap- 
proach is completely justifiable. On the other hand, the 
same process can be addressed in the target rest frame 
[Tol [Tl| . In this framework, the direct photon produc- 
tion is viewed as an electromagnetic bremsstrahlung of 
a (light or heavy) quark that interacts with the target 
via gluonic exchanges. This idea has been used to de- 
scribe Tevatron data on prompt photon production at 
midrapidities using light quarks contribution 0, 0] ■ In- 
vestigations of such formalism in association with aspects 
of saturation physics, for instance the Color Glass Con- 



densate framework fT2| , have been made for dA and AA 
collisions in a number of works p"3l - [T5j . 

In this work we apply the target rest frame formalism 
to photon plus heavy quark production cross section in 
high energies. The dependence of the production cross 
section on the final state particle rapidities and trans- 
verse momenta is investigated. In next section we intro- 
duce the main formulae to compute the hadronic differ- 
ential cross section. It is already known from previous 
studies Q that the photon transverse momentum spec- 
trum is a sensitive probe of color dipole scattering cross 
section. We investigate the consequence of using several 
implementations of dipole cross sections which are con- 
strained from deep inelastic scattering data. The theo- 
retical results are compared to Tevatron measurements 
for the 7 -|- Q-\- X production cross sections. An analysis 
for proton-proton and proton-lead collisions at the LHC 
regime is also developed. In order to do so, we will rely 
on geometric scaling arguments to do the transition from 
a nucleon to a nucleus target. The numerical results are 
discussed in detail in last section. 



II. PHOTON IN ASSOCIATION WITH HEAVY 
QUARKS IN HIGH ENERGY HADRON 
COLLISIONS 



In this section we summarize the relevant formulae to 
compute the Q -I- 7 production cross section at small- a; 
regime in the target rest frame. Such an approach has 
been early developed by Kopeliovich [Ifl] to the Drell-Yan 
production in the context of nuclear shadowing and fur- 
ther developed in [lH . Here, it will be considered that the 
photon transverse momentum is large compared to the 
hadronic scale, Aqcd, and at the same time smaller than 
the hadronic center-of-mass energy, Ag^^ ^ "p^^ <^ s. 
In the target rest frame the production process is viewed 
as follows. A large- a; heavy quark (or heavy antiquark) 
of the hadron projectile scatters off the gluonic field of 
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the hadron target and radiates a real photon. The rel- 
evant diagrams are those where the photon is radiated 
before or after the interaction with the target, whereas 
diagrams involving the quarks interacting with the target 
both before and after the photon vertex are suppressed 
in the high energy limit. Although in the process of elec- 
tromagnetic bremsstrahlung by a heavy-quark no quark- 
antiquark dipole participates, we soon observe that the 
cross section can be expressed via the elementary dipole 
cross section, Ngg{x,r,b), of interaction of a qq dipole 
with a target. This allows us to make use of the exten- 
sive phenomenology on deep-inelastic scattering where 
the dipole cross section is very well determined from data 
on small- a; region [l6l |. 

The scattering cross section for production of a mass- 
less on-shell quark with momentum I and a real pho- 
ton with momentum k was derived in [H M- The 
differential cross section for the process Q{p) + h — > 
[Q{l)j{k)] +X is given by: 



is related to the photon and final state quark rapidities 
and transverse momenta as follows ^14i |: 
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In order to compute the hadronic cross section for pho- 
ton plus heavy quark production it is needed to convo- 
lute the partonic cross section, Eq. ((TJ, with the heavy 
quark distribution function (PDF) on a proton, fq^x, Q). 
Therefore, the differential cross section for the hadronic 
production is given by. 



d^a {pp ^ 7Q -t- X) 
d{p\Y dy^ dyq 
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d{p\Y dy^ dyq' 
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with the quantity H being defined as, 
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where the incoming heavy quark has momentum p, the 
photon and outgoing quark rapidities are defined via 

= fc- = ^e^^ and /3q = r = ^e^^Q. The cor- 
responding transverse momenta of photon and outgoing 
quark are p^ and p^, respectively. The angle At/i is the 
opening angle between the final state quark and photon 
defined as cos (A(/)) = {jyC ■ PKl^) /p'Yp^, with respect 
to the produced quark axis. The quantity N is related 
to the color dipole cross section, N{x,r), which should 
satisfy the QCD evolution equations and that includes 
the small- a; evolution. It is defined as: 



AfiXg.pV^pl'^ 



d^b / d^re'^-P"^N{xg,r,b), (3) 



(7ri?2) 2T:rdrJo{rpT)N{xg,r), 



The charm and bottom PDFs are presently assumed to 
be radiatively generated and then are related to the gluon 
PDF through DGLAP evolution equations. It is noticed 
that the charm structure function at large-x measured 
by EMC Collaboration suggests that there might be an 
intrinsic charm (IC) component in nucleon. In general, 
nonperturbative models are considered to describe such 
(l)a component (see, for instance Refs. fT7l-(T9j). Those 
models are not considered here. In Eq. ([S]), the differ- 
ence between charm and bottom comes from the quark 
charge, = 4/9 and el = 1/9, and from the quark PDF 
(charm PDF is larger than the bottom one). For the 
factorization/renormalization scale, Q — fi, the photon 
transverse momentum is usually chosen. 

To proceed further, we need to know the quantity 
J\f{xg,p±'*,pj_'^). It can be computed numerically us- 
ing the available phenomenological dipole cross sections 
|16| . Here, we consider some of them that have been 
used to describe deep inelastic scattering (DIS) data and 
RHIC data as well. We can gain some physical insight 
using analytical versions of the amplitude Af. For in- 
stance, considering small dipole configurations (which is 
the typical configuration for large px considered in the 
current study) the color dipole amplitude behaves like 
N{r — > 0) « (r'^Qs)'' /4:. Such an expression is typical 
in saturation models based on the QCD nonlinear evolu- 
tion equations for the dipole amplitude and its behavior 
is known as geometric scaling on the variable rQ^. The 
intrinsic momentum scale, Qs oc x^^, is the so-called sat- 
uration scale (For central rapidities at Tevatron, Qs ^ 1 
GeV). For some fixed values of the anomalous dimension, 
7, the zeroth-order Hankel transform in Eq. can be 
done analytically considering pr > 0- In Ref. |14| . it has 
been shown that for 7=1/2 (BFKL anomalous dimen- 
sion) the result is the following: 



where the b is the impact paramater of interaction and r 
is the corresponding transverse size. Here, we have used 
the notation = (pj_^+P±'^) andpT — \pt\- In last line 
in Eq. ^ we treat the target (proton) as a homogeneous 
disk of radius Rh — 5 GeV. The momentum fraction Xq 



U{xg,pV,p\'^)_ ^ 327r 
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The current Tevatron data on 7 -|- Q production is dom- 
inated by large photon transverse momentum, p\ > 30 
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GeV ^ Qs, and then it is theoretically expected that 
the anomalous dimension is close to the DGLAP values, 
7DGLAP = 1- Thus, it is important to consider dipole 
models where the anomalous dimension is running. This 
is named extended geometric scaling property. For a run- 
ning anomalous dimension which does not depend explic- 
itly on the dipole size r, it can be obtained the following 
[20]: 
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where we consider two models for the running anomalous 
dimension, which have been confronted to RHIC data on 
charged particle multiplicities. The first one is the BUW 
model ^2^, whereas the second one is the DHJ model 
plj . The expression for 7(0;), with uj = px/Qs, for each 
model is given by: 



7BUW 



7dhj 



7. + (1-7.) 



7. + (1 - 7s) ■ 
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The BUW parameters are a = 2.82 and 6 = 168 [20| 
and the DHJ one is d = 1.2 (in addition, one has Y = 
ln(l/a;)) [2T|. In both cases, 7^ = 0.628, A = 0.3 and 
xq = 3 • 10"'^, with a saturation scale on the nucleon 
given by ^ (xo/xg)^ GeV^. 

Furthermore, in Ref. [l^ it has been shown that the 
quantity J\f is directly connected to the unintegrated 
gluon distribution, J-(x,k±)^ in the small dipole size 
(r — >■ 0) limit. Using such a connection, we can now 
use the available numerical/analytical implementations 
of the unintegrated gluon function. In next section, we 
calculate the differential cross section, Eq. ([5]) , using the 
analytical expressions for the dipole amplitude, Eq. ([T]), 
computed from BUW and DHJ models and a numerical 
result for the unintegrated gluon distribution. We will 
discuss the main features coming from distinct theoreti- 
cal approaches when compared to the NLO perturbative 
QCD calculations. A study for an extrapolation to pp 
and pA collisions ate the LHC will be also presented. 



III. DISCUSSIONS AND SUMMARY 

Let us now compare the theoretical calculations pre- 
sented in previous section to the experimental data for 
production cross sections in pp collisions at Tevatron for 
center of mass energy ^/s = 1.96 TeV. Recently, DO 
Collaboration have performed the first measurement of 
the differential cross section of inclusive photon produc- 
tion in association with heavy flavors jets [s^. The re- 
sults cover the range in photon transverse momentum 
30 < p\ < 150 GeV ad photon/jet rapidities, \yj\ < 1 
and |j/jot| < 0.8. The following cut is required to jet 
transverse momentum > 15 GeV. For simplicity, here 
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FIG. 1: The 7 + c + X differential cross section as a function 
of p'Y- Experimental data from DO Collaboration 8] (only 
systematic errors are presented). 
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FIG. 2: The j + b + X differential cross section as a function 
of p\. Experimental data from DO Collaboration 8j (only 
systematic errors are presented). 



we will take central rapidities y^, yq — and consider, 
2/jct — yq- Concerning the numerical calculation of dif- 
ferential cross sections in Eq. ([IJ some comments are in 
order. As already discussed in Ref. [l3|, despite any par- 
ticular assumptions on the dipole amplitude the factor 
in denominator of Eq. ^ diverges as the momenta of 
produced quark and photon are parallel, A0 = 0, which 
is the usual coUinear divergence present in perturbation 
theory. Thus, in our calculations we consider a lower 
cut-off for the angular integration, A</)iiiin- Moreover, a 
numerical calculation of the Hankel transform of zero or- 
der in Eq. ^ is divergent when \p±' = 0, which 
occurs for A0 = tt if momenta are equal, p^ = p^. In our 
case, this is not a shortcoming as analytical expressions 
for the color dipole amplitudes are being used. For a 
similar discussion connected to the present study on the 
angular dependence of production cross section we quote 
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™" < 0.08 



the recent work in Ref. [2j|. As a last comment, we are 
using the target rest frame approach in the hmit of its 
vahdity in the Tevatron regime. Formally, this approach 
is valid at small- a;^, see Eq. and a careful analysis 
on what is the a:g-range being probed is difficult as the 
quark transverse momentum, p^, is integrated over in 
Eq. (HI). However, the minimum Xg value can be esti- 
mated. Using the Tevatron quark transverse momentum 
cut, > 15 GeV and the photon transverse momentum 
range, 30 < pX < 150 GeV one gets 0.02 < x 
at central rapidities. 

In Figs. [T] and [5] the 7 + c + X and 7 + 6 + X dif- 
ferential cross section arc presented as a function of p''^ . 
Experimental data from DO Collaboration ^8,] are shown, 
where only the systematic errors are presented. The long- 
dashed lines (in green) represent the the NLO pQCD pre- 
dictions using CTEQ6.6M PDFs [H] (averaged over 
rapidity regions), where the renormalization scale ^ij, 
factorization scale fip and fragmentation scale ^/ are all 
set equal to p\. In our numerical calculations the same 
set is used, i.e. in Eq. ([5]) we consider /^^ = {p\Y ■ 
The long-dashed curves (in blue) correspond to the cal- 
culation using the connection of the dipole amplitude M 
to the unintegrated gluon function. For the unintegrated 
gluon PDF we consider the J. Bluemlein parametrization 
(JB) [251], which relies on BFKL approach [111: 

J Xg Tj Tj 

where it is defined the Q function as, 



Qiri.vl) ^^Jo f2.\/a,ln(l/ry)ln(MVp|)) , (Pt < m') 
VPt \ ^ J 

G{v,Pt) ^^Io (2JasH^/v)HPT/^^^)] : (pt > m'), 
VPt \ / 

where Jq and Iq stand for Bessel functions (of real and 
imaginary arguments, respectively) , and as — as/i-K. 
The LO MRST set [23 was used in our calculations as 
the input gluon coUinear density, g{x, Q). 

Let us now comment on the numerical results for the- 
oretical approaches considered here. As already found in 
Ref. 0, the NLO pQCD calculation describes in good 
agreement the bottom data. However, the large range 
of charm data is underestimated and it seems to favor 
a IC component which would be important at that re- 
gion. On the other hand, the target rest frame approach 
produces distinct results depending on the input for the 
dipole cross section. The JB parameterization is close 
to the NLO pQCD (the anomalous dimension for JB is 
near the DGLAP one as it is derived in the double log- 
arithmic, DLL, limit of BFKL approach) for charm and 
deviates from pQCD at p^ > 60 GeV for the bottom case. 
A weak reason could be due the fact we are considering 
the strong coupling fixed, = 0.2, in JB parameteriza- 
tion for both cases and a distinct bottom PDF. We have 
checked that it is not the case. Probably, the deviation 
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FIG. 3: The j + Q + X differential cross section at LHC 
regime. Left panel: proton-proton cross sections at central 
jet /photon rapidities as a function of at = 14 TeV. 
Right panel: proton-lead (pPb) cross section for charm plus 
photon for fixed photon rapidities y~' = 1, 2, 3, 4 at ^/s = 8.8 
TeV. 



does come from the relative contribution of a set of dia- 
grams that are not accounted for in the target rest frame 
approach compared to pQCQ. For instance, it is shown in 
Ref. [§| that the annihilation process qq — QQ"f is very 
important at large p'[_ for the bottom case. The target 
rest frame approach only involves the Compton subpro- 
cess and can not be directly compared to the NLO pQCD 
calculation. 

The BUW and DHJ parameterizations differ not only 
in the scaling behavior, but in the way the large px limit 
of 7(0;) approaches to one. This last feature leads to dif- 
ferent large momentum slopes of dipole amplitude and 
therefore to distinct predictions for the large pT -slope 
when using Eqs. ([S]), ^ and ([7]). For instance, in [20| 
it was shown that for BUW model, A/buw(7 1) oc 
[Qs/ptY^"' /Pt^ whereas for DHJ one has A/dhj(7 



1) cx Qr°/bTln(p|/Q|)]. Therefore, the future data 
on forward production (where current approach is for- 
mally valid even at large pY) of heavy quarks in associ- 
ation with prompt photon could discriminate among the 
distinct models for the dipole cross section. 

Finally, we perform some estimates for the LHC kine- 
matic regime. Let us start by pp collisions at energy 
of = 14 TeV. Now, even at central rapidities the 
Xg is sufficiently small to justify the target rest frame 
approach at large transverse momenta. Keeping the 
same kinematic cuts from Tevatron, > 15 GeV and 
30 < < 150 GeV, one gets 0.003 < a;™ < 0.01. 
The situation is further improved in the forward rapidi- 
ties case. In Fig. [3] (left panel) are shown the estimates 
for the differential cross section (for charm and bottom) 
as a function of photon transverse momentum at cen- 
tral rapidities. The dashed curves represents the results 
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for BUW model, whereas the dot-dashed ones stand for 
the DHJ model. The deviations at large p]_ follow the 
same pattern as for Tevatron. The conclusions are similar 
as previously discussed, where the transverse momentum 
distribution allows discrimination among phenomenolog- 
ical/theoretical models for the dipole scattering cross sec- 
tion (or for the unintegrated gluon distribution in an in- 
direct way). At the LHC the center-of-mass energy is 
sufficiently high and annihilation process no longer dom- 
inates at large p\ as seen for Tevatron, mostly for the 
bottom case. 

In Fig. [3] (right panel) we present the proton-lead 
(pPb) differential cross section at ^/s = 8.8 TeV as a 
function of photon transverse momentum for fixed for- 
ward photon rapidities y'^ = 1,2,3,4 (from the top to 
bottom, respectively). In order to obtain the cross sec- 
tion for a nuclear target we rely on the geometric scaling 
arguments [28|: it is replaced Rp — ?> Ra in Eq. ([3]) and 
also Ql p {ARl/R\)^ Ql ^, where A = 1.26. In case 
of A = 1 that replacement becomes the usual assumption 
for the nuclear saturation scale, Q?4 = ^^^^ Qs,p- Such 
an approach has been used in Ref. [20| to describe small- 
X data on the nuclear structure functions. This enhances 
the saturation scale by a factor six for a lead nucleus. 
As photon rapidity diminishes, i.e. smaller Xg, the satu- 
ration scale increases since Qs oc x~'^. Such a behavior 
implies an enhancement on the quantity w = pt/Qs at 
fixed p]_ in Eq. ([7]) and therefore large rapidities means 
bigger w and the running anomalous dimension increases 
towards 7 — >■ 1. This patters is viewed in the plot corre- 
sponding to a steep p]^-slope as rapidities increase. We 
have checked that the p]^-slope is unchanged by consider- 
ing forward quark rapidities, t/jet = yg ^ 1, at fixed cen- 
tral photon rapidity. As a final comment, the 7 -f Q Pro- 
duction in pPb is important per se as the nuclear gluon 
distribution is largely unconstrained and such a process 



can test at the same time the gluon and charm/bottom 
distribution functions. Thus, measurements having ap- 
propriated experimental error could distinguish the dis- 
tinct parameterizations for the nuclear PDFs. 

As a summary, we have applied the target rest frame 
formalism to photon plus heavy quark production cross 
section in hadronic collisions at high energies. The de- 
pendence of the production cross section on the photon 
and quark rapidities and transverse momenta have been 
investigated. It was verified that the photon transverse 
momentum slope is a good probe of color dipole scat- 
tering amplitude as it directly depends on the (running) 
anomalous dimension, which is determined from the un- 
derlying QCD dynamics in a given kinematic regime. 
The theoretical results are first compared to Tevatron 
measurements of the differential 7 -I- Q+ X production 
cross sections having in mind that this is an extrapola- 
tion on kinemetic regime where the small-Xg limit id not 
completely fulfilled. An analysis for proton-proton and 
proton-lead collisions at the LHC regime, including the 
forward rapidity region, has been performed. For the 
pA model for the saturation scale for a nucleus 

has been introduced based on geometric scaling argu- 
ments, where Ql{xg,A) oc A^^^Q^Xg). This procedure 
can be also done using a nuclear version of the uninte- 
grated gluon distribution as proposed for instance in {29{ . 
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